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In this paper we present generic properties of quantum transport in mono-layer graphenc. 
In the scheme of the Kubo-Geenwood formula, we compute the square spreading of wave 
packets of a given energy with is directly related to conductivity. As a first result, we 
compute analytically the time dependent diffusion for pure graphene. In addition to the 
semi-classical term a second term exists that is due to matrix elements of the velocity 
operator between electron and hole bands. This term is related to velocity fluctuations i.e. 
Zitterbewegung effect. Secondly, we study numerically the quantum diffusion in graphene 
with simple vacancies and pair of neighboring vacancies (divacancies), that simulate 
schematically oxidation, hydrogenation and other functionalisations of graphene. We 
analyze in particular the time dependence of the diffusion and its dependence on energy 
in relation with the electronic structure. We compute also the mean free path and the 
semi-classical value of the conductivity as a function of energy in the limit of small 
concentration of defects. 
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1. Introduction 

Graphene consists of a mono-layer carbon atoms forming a 2D honeycomb lattice. 
It has attracted much inter est b ecause of its unique electronic properties and poten- 
tial for device applications Because of the linear dispersion relation of electron 
states close to Fermi energy, many unconventional transport behaviors have been 
shown in graphene. The effects of defects and impurities are essential in trans- 
port properties. In particular the role of d efect a nd impurities in a Metal-Insulator 
transition has been found experimentallypEEl but is still not completely under- 
stood theoretically. 4 ^ 13 Atomic vacancy (Va) is one of the natural d efects pre sent 
in graphene that has been intensively studied theoretically (see Ref. ' ' ' ' ' and 
Refs. therein). Moreover, as shown recently by ab initio studies pi vacancies simu- 
late schematically oxidation and chemical adsorption of simple atoms or molecules 
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(H, OH, CH3...). Indeed, adsorption of atom creates a covalent bond between an 
orbital of this atom and a p z orbital of one carbon of the graphene sheet. There- 
fore, this p z orbital does not contribute any more to the conduction band (p z 
band). Thus, roughly speaking, the effect on conduction of one adsorption band 
is simulated by removing one C atom from graphene sheet. In the same scheme, 
an adsorption that c reat es two covalent bounds on two neighboring C atoms (for 
instance epoxy group ^^J) are schematically simulated by a pair of two neighboring 
vacancies ( "vacancy pairs" or divacancies) . 

In this paper, we present theoretical investigations of transport properties in 
graphene. In the second part, we discuss the link between conductivity and quantum 
diffusion. In the third part, it is shown that inter-band transitions have significant 
effect on quantum diffusion in pure graphene. In the fourth part, effects of simple 
vacancies and divacancies on quantum diffusion are analyzed. 



2. Conductivity and quantum diffusion 

In the framework of Kubo-Greenwood approach for calculation of the conductivity, 
a central quantity is the average quadratic spreading of wave packets of energy E 
at time t along the x direction, 



AX 2 (E > t) = ( (X(t)-X(0)) ) , (1) 



where X (t) it the Heisenberg representation of the position operator X. (A) e means 
an average of diagonal elements of the operator A over all states with energy E. 
The diffusivity at zero temperature, T>(E), at energy E is deduced from AX 2 , 

V(E)= lim D(E, t) with D(E,t) = AX2( - E ' f \ ( 2 ) 

t— >+oo t 

where D(E, t) is called diffusion coefficient. In a 2-dimensionnal system with surface 
S, the DC-conductivity a xx at zero temperature along the ^-direction is given by 
Einstein formula: 

e 2 

i xx (Ef) = -n(E F )V(E F ), (3) 

where n(E) is the total density of states and E-p the Fermi energy. 

The effect of decoherence mechanisms such as electron-electron scattering, elec- 
tron phonon interaction (temperature), is not considered in the above expression. 
This effect can be estimated by introducing an inelastic scattering time r, ; . Tj de- 
creases when the temperature incr ease s. In actual graphene at room temperature, 
realistic values of Tj is a few 10 -13 sP^ Therefore the experimental conductivity can 
be estimated by: 

e 2 

<J xx (Ef) - -n(E F )V(E F , Ti ), (4) 
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The variation of the experimental conductivity with t\ is analogous to the variation 
with time which is discussed in detail in this article. Here the Fermi-Dirac distri- 
bution function is taken equal to its zero temperature value. This is valid provided 
that the electronic properties vary smoothly on the thermal energy scale k B T. 



3. Quantum diffusion in pure graphene 

We have develope a simple tight binding (TB) scheme that reproduces the ab 
initio calculations of the electronic states with energies ±2eV around the energy of 
Dirac point {E = 0). Only p z orbitals are taken into account since we are interested 
in what happens at the Fermi level. The Hamiltonian is, 



(5) 



<hj> 



where coupling matrix element, = (i\H\j), depends on the distance 
orbitals \i) and \j), 



Tij between 



Uj = - 7 oe^ (1 - r - /a) , 



(6) 



with a = 1.418 A, the firt neighbor distance. First neighbors interaction in a plane 
is taken equal to 70 = 2.7 eV. Second neighbors interaction 7q in a plane 

equalsPl 

to 0.1 x 70 which fixes value of the in (j6|). 

In crystals, once the band structure is calculated from the tight-binding Hamil- 
tonian the a ver age quadratic spreading can be computed exactly in the basis of 
Bloch states ff^ In pure graphene the average square spreading is the sum of two 
terms ff^ 



AX 2 (E,t) = Vit 2 



AX 2 B (E,t). 



(7) 



The first term is the ballistic (intra-band) contribution at the energy E. Vb is the 
Boltzmann velocity in the x direction. The semi-classical theory is equivalent to 
taking into account only this first term. The second term (inter-band contributions), 
AX^ B (E, t) , is a non- ballistic (non-Boltzmann) contribution. It is due to the non- 
diagonal elements in the eigenstates basis {\n)} of the velocity operator V x , 



ih 



X,H 



From the definition (TTJ), one obtains 



H5I2DI 



(8) 



AX 2 B (E,t) = 2h 2 



E 



* k ,n' '(n J '^n) - 



1 — COS 



(nk\V x \n'k) 



Er =E 



where E^ n is energies of eigenstate \nk). 
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In graphene, assuming a restriction of the Hamiltonian ([5]) to the first neighbor 
interactions only, AX^ B (E,t) is given by: 

AX^(i^) = ^(l-cos^) withy B = ^|. (10) 

At small time, t < h/E, the Non-Boltzmann term is equal to Boltzmann term 
AX£ B (£,i) ~ V 2 t 2 , thus AX 2 (E,t) ~ 2V; 2 i 2 and £)(£,*) ~ 2V B 2 i . Whereas for 
large t, the Boltzmann term dominates and AX 2 (E,t) ~ V^i 2 and D(E,t) ~ V^t. 
The non-Boltzmann term is due to matrix elements of the velocity operator be- 
tween the two bands (i.e. between the hole and electron states having the same 
wavevector). These matrix elements imply that the velocity correlation function 
has also two parts: one constant and the other oscillating at a frequency 2E/H 
where E is the energy of the state. This is precisely the phenomenon of jittery 
motion also called Zitterbewegung. Note that in any crystal having several bands 
there are also components of the velocity correlation function which are oscillat- 
ing at frequencies {E^ n — E^ n ,)/fr- Therefore Zitterbewegung is quite common in 
condensed matter physics. For example approximants of quasicrystals present very 
strong Zitterbewegung effe ct and the non-Boltzmann contribution dominates the 
Boltzmann contributionP222Il22] 

4. Quantum diffusion in the presence of vacancies 

In this part, the density of states (DOS) and diffusion coefficient calculated nu- 
merically in disorder graphene are presented. We study both the effects of sim- 
ple vacancies and the effects of pairs of nearest neighbor vacancies (divacan- 
cies). The total DOS is computed by a Lanczos-type recursion method. In the 
scheme of the Kubo-Greenwood formula for transport properties, the diffusion 
coefficient, D(E,t), is computed by using the PEQD method (Polynomial Ex- 
pansion for Quan tum Diffusi on method) developped by the Mayou, Khanna, 
Roche and Triozoni^SEHZSET] This method allows very efficient numerical calcula- 
tions by recursion in real-space. This method has been used to study quantum 
transport in disordered graphene^ and chemically doped graphene graphene 
with Ozone functionalizationpl graphene with structural defects (Stone- Wales 
and divacanciesj^. We use periodic boundary conditions with a supercell size 
L x = 2947.2 nm, L y = 122.8 nm. The energy resolution is about 0.01 eV. Defects 
(simple vacancies and divacancies) are randomly distributed in a supercell. We have 
checked that the present results do not depend on the size of the super cell. 

Electronic structure of simple vacancies has been intensively studied JHEIEEIIUZI 
Graphene lattice is a bipartite lattice with two equivalent sublattices A and B. If we 
consider an Hamiltonian containing only nearest-neighbor hopping (uncompensated 
latticePS) 

each A atom is coupled with only B atoms (and reciprocally). Conse- 
quently, at very low concentration of vacancies (~ isolated simple vacancies), the 
distribution of vacancies is loc ally u neven between the two sublattices and zero en- 
ergy states necessarily appears I^SIZI But if vacancies are arranged by pair of nearest 
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Fig. 1. DOS in graphene for various concentration of simple vacancies and divacancies. 
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Fig. 2. Diffusion coefficient, D = AX 2 /*, versus time t for wave packets at energy E in the peak 
of simple vacancies: E = — 0.23eV. 



neighbor vacancies (divacancies), the electronic structure at low impurities concen- 
tration is completely different. Indeed in that case, the distribution of vacancies is 
locally even between the two sublattices, and zero energy state does not occur. 



I 



I 
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Fig. 3. Diffusion coefficient, D = AX 2 /t, versus time t for wave packets at energy E, with various 
concentrations of vacancies: (point line) perfect graphene, (full line) divacancies, (dashed line) 
simple vacancies. 



A finite concen tratio n c of simple vacancies induced a small enlargement of the 
zero energy modesPZH Moreover, for Hamiltonian with second neighbor hop ping, 
these states are strongly enlarged and their energy is shifted to negative values! I I 
As shown figure [1] our TB model reproduces this behavior in good agreement 
with previous studies. The states due to simple vacancies is observed at the energy 
E p ~ — 0.23eV. As expected, the peak at E v ~ — 0.23eV is not obtained for small 
concentration of divacancies (figure [lj . For larger concentration of divacancies the 
distribution of vacancies could be locally uneven between the two sublattices and 
the peak at E p ~ — 0.23eV is found. For large concentration of simple vacancies 
and divacancies, c > 30%, a low density of states is obtained around E = and 
gaps may occur for c > 50%. 

The diffusive coefficient D, D = AX 2 /t, versus time t, is shown on figure [5] 
for an energy in the peak of simple vacancies, E — —0.23. Different regimes are 
observed. For small time t, the effect of elastic scattering is very small and the 
regime is ballistic, D(E,t) = V 2 t. As expla in in the previous section, the effective 
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Fig. 4. Semi-classical conductivity a sc (in unit of Go = 2e 2 /h) and elastic mean free path l e versus 
energy E: (a) and (c) simple vacancies, (b) and (d) divacancies. Concentration of defects: cercle 
c = 0.01%, triangle c = 0.1%, cross c = 1%, diamond c = 10%. 



velocity is V = V2Vb for t < h/E and V = Vb for t > h/E, which leads, in this 
log-log scale, to a shift of D(t) at t c± h/E (figures [5] and [3]). After the ballistic 
regime D reaches a maximum values equals to the semi-classical diffusivity T> SC (E). 
Indeed if there were no localization effects the diffusivity would be constant at 
larger times. The elastic scattering time r e {E) is defined by r e = T> sc /Vq, and the 
elastic mean free path, l e , is deduced from l e = Vqt c . At large time, in 2D-system 
with defects, quant um interferences lead to decay of D(t) with a characteristic 
localization length £p22 When l e -C the semi-classical diffusive regime (saturation 
of D, D(E,t) ~ D SC (E)) occurs on a large time scale. But for l e ~ £, the semi- 
classical approximation breaks down and the decay of D(t) is strong for t > r e . 

The mean free path and semi-classical value of the conductivity are shown in 
figure S] as a function of energy and for different concentrations. The semi-classical 
value of the conductivity is a sc (Ep) = ^■n(Ep)T) sc (Ep), where T> sc (Ep) is semi- 
classical diffusivity defined above. Let us recall that when the semi-classical of 
conductivity is much larger than the quantum of conductance, Gq — 2e 2 /h, the 
localization length is much larger than the mean free path and there is a large scale 
of inelastic scattering times Tj on which the quantum interferences are small and 
the semi-classical theory of transport is valid. 
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Fig. 5. Diffusion coefficient, D = AX 2 /t, versus energy E of wave packets at time t = 3.26 X 10 s, 
for (a) c = 0.1% vacancies, (b) c = 10% vacancies. 

One sees from figure|2]that for sufficiently low doping the conductivity and mean 
free path are inversely proportional to the concentration of defects in agreement 
with the Fermi Golden Rule. For divacancies the mean free path increases when 
the energy tends to zero. This is in agreement with the semi-classical model since 
it can be shwon that the T-matrix of the divacancy is energy independent close to 
the Dirac point and therefore th mean free path is inversely proportional to the 
energy. For simple vacancy, it is also clear that close to the energy of the resonant 
state E = — 0.23cV the scattering is strong and accordingly the semi-classical 
conductivity and mean free path are small. In the energy region —0.4 < E < 0.4 eV 
which is of most interest for experiments the mean free path can vary by a factor of 
about ten. This variation is partly compensated by the variation of the density of 
states so that the semi-classical conductivity varies by a factor of about 3 at most 
in the same energy range. 

The diffusive coefficient D, D = AX 2 /t, versus the energy E of the wave packets 
is shown on figure [5] for two concentrations of defects. Figures |3] and [3] compare 
diffusion coefficient D in both case: simple vacancies and divacancies. For energies 
out of the peak in the DOS due to simple vacancies (E = and E = 0.23 eV in 
figure [3]), this two kinds of defects have similar consequence on diffusion. But for 
energy in the peak [E ~ — 0.23eV) and defect concentration c < 10%, electron 
localization by simple vacancies is stronger. This is due to resonant scattering by 
the impurity potential. This is very similar to scattering by virtual bound states 
in the physics of transition metals for example. As shown figure [5jb), for large 
concentration of defects (c > 10%), simple vacancies and divacancies have similar 
scattering effects. 

5. Conclusion 

To conclude we have studied in detail the quantum diffusion of electrons in graphcnc 
in the presence of local defects. These local defects are simple vacancies and diva- 
cancies. They are representative of functionnalizcd graphcnc with molecule that 
produce covalent bonding either with one carbon atom (simple vacancies represent 
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for example fonctionalization by Hydrogen)) or with two neighboring carbon atoms 
(divacancies can represent functionalization by one oxygen atom (epoxy group)). 
The time dependence reveals different regimes: ballistic with a kink at short time, 
saturation, existence of a plateau if the localization length is much greater than the 
mean free path, decrease of the diffusivity due to localization. We show also impor- 
tant differencies between the electronic structure of simple vacancies and pairs of 
nearest neighbor vacancies (divacancies) . In the case of simple vacancies a narrow 
maximum occurs in the density of states due to the existence of an impurity res- 
onance. Just at the impurity energy (E ~ — 0.23eV) the scattering is strong and 
the tendency to localization is enhanced. Finally, we present also the mean free 
path and the semi-classical value of the conductivity as a function of energy in the 
limit of small concentration of defects. Let us recall that when the semi-classical 
value of the conductivity is much larger than the quantum of conductance 2e 2 /h 
the localization length is much larger than the mean free path and there is a large 
scale of inelastic scattering times Tj on which the quantum interferences are small 
and the semi-classical theory of transport is valid. 
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